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It is hard to avoid the scanning tunneling microscope 
(STM) in the imaging, fabrication, or characterisation of 
nanoscale structures today . It is certainly the most 
ubiquitous tool in surface science. But while the for- 
mation of atomic wires from metal leads and their final 
breaking is well researched j^Q , the physics of the reverse 
process - an STM tip approaching the surface at very 
close range - is far less well documented. This is so even 
though the phenomena that occur - the forces, atomic re- 
laxations, the sudden changes in image at a critical point 
- are the everyday experience of every STM experimenter. 
Even more importantly, several routes are currently be- 
ing explored to fabricate and operate electronic devices 
on the nanoscale . In chemical methods of nanofab- 
rication || a high precision of atomic assembly is attain- 
able, even though the exact position of every atom is not 
controlled. By contrast, to initiate such a process, or to 
craft structures by moving each element into place, the 
exact position of the STM tip and the target atoms must 
be known. Furthermore, nanoscale research requires a 
high degree of knowledge about the electronic, chemical, 
and transport properties of potential constituents - in- 
formation often obtained by STM. This information may 
be misleading if the tip-sample distance is not correctly 
estimated. Therefore, as soon as the work moves from 
qualitative to quantitative, distance estimates are of crit- 
ical importance. The question of distances has previously 
been analysed to some extent for STM &- p||, and for the 
related scanning force microscope (SFM) |ig-|l|. What 
is missing, so far, is a picture of the interplay between 
forces, relaxations, and currents. 

In this paper we analyze the interactions between tip 
and sample and their effect on distance estimates and 
STM images from first principles. We demonstrate that 
our simulation technique allows us to keep track of all 



key quantities in this situation. We choose the Au(lll) 
surface since it is widely used in nanotechnology and was 
studied in recent extensive experiments [jl3 14 1. 

The computational method is based on first principles 
density functional theory. As is standard practice in cal- 
culations of this type only a small part of the periodic 
surface is used to determine the physical properties of 
the whole system. Moreover we employ periodic bound- 
ary conditions both parallel and perpendicular to the sur- 
face. The sample surface was mimicked by a five layer 
Au(lll) film with nine atoms per layer. The high number 
of atoms per layer was necessary for a complete decou- 
pling of the interactions between adjacent tip pyramids. 
The central layer of the Au(lll) film was our plane of ref- 
erence. The tip apex was mimicked by a tungsten tetra- 
hedron mounted on the reverse of the Au(lll) surface 
slab; the tip-sample separation was altered by varying 
the size of our unit cell normal to the surface. Tungsten 
is the most commonly used tip material in experiments; 
the low reactivity and the similar size of tungsten and 
gold atoms makes the support of the tip apex on gold 
realistic. We also checked its plausibility by a careful 
analysis of the forces during our simulations. In this case 
the force gradients are mainly due to the position of an 
atom, but not its chemical nature. In Fig. [l] we display 
the forces on the apex atom and the subsurface atoms of 
the tip, it can be seen that the forces on the apex atom 
are about one order of magnitude higher. 

Even for the pseudopotential VASP (Vienna ab initio 
Simulation Program Jl^,|l6|) code the simulations were 
rather demanding and at the limit of realistic time scales. 
This fact ruled out a check of the results including addi- 
tional gold or tungsten layers. For the same reason we 
did not include long-range forces. But as established by 
Perez et al. in their analysis of atomic force microscopy 
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PI , they only cause significant relaxations in the very 
low distance regime (below 5 A), and in this range short- 
range forces are substantially higher, as we find in our 
simulations. The same point was made in an earlier pub- 
lication by Diirig et al. | 17|] . We also neglect bias-induced 
surface relaxations in our treatment. The technical de- 
tails of the calculation are the following: The cutoff of 
the ultrasoft pseudopotentials JlJ] was 230 eV, exchange 
and correlation effects where computed with generalized 
gradiant correction and using the formulation of Perdew 
et al. |[{|. The Brillouin zone was sampled with (4x4 
x 1) k-points of a Monkhorst-Pack grid p0|, forces were 
converged to values of less than 0.01 eV/AT 

The reading of a piezoscale in the experiments is given 
by the position of the plane of reference in the supercell 
setup of our system. From an initial tip-sample separa- 
tion of 10 A the distance was gradually reduced, initially 
in steps of 0.5 A. At every step the system was fully re- 
laxed. The onset of significant mechanical interactions 
between tip and sample occurs at about 5.0 A. To de- 
termine the precise effect of these interactions in the low 
distance range we reduced the distance in steps of 0.1 A 
from 5.0 to 4.0 A for one tip position (STM tip on top 
of a surface atom), and from 4.5 to 3.5 A for the second 
tip position (STM tip between the surface atoms). We 
computed the surface wavefunctions for the fully relaxed 
sample and tip systems at every distance and for two tip 
positions. In this part of the calculation the initial sys- 
tem was separated into two subsystem, the sample and 
the tip surface, and the surface wavefunctions computed 
for each subsystem separately. These wavefunctions were 
used as the input for our simulation of the tunnel currents 
pl[ . The procedure leads to separate constant current 
contours of the surface, depending on the position of the 
tip. Interpolating between the two simulations we obtain 
the change of corrugation due to system relaxations from 
first principles. 

It has been shown by a number of authors [^,| 1 1 1 that 
substantial relaxations are confined to the low distance 
regime. Perez et al. obtained their result with an ab ini- 
tio method, but they calculated the changes on a semi- 
conductor surface and did not relate them to measured 
currents. Olesen et al., and Clarke et al. analysed the 
events on metal surfaces, but employed empirical poten- 
tials; for this reason they too could not simultaneously 
evaluate the current. Their results show, not surpris- 
ingly, that the onset of forces and relaxations begins 
when the interatomic distances are close to the typical 
bondlength of the atoms. However, from experiments 
we know that tip-sample forces on Au(lll) are apprecia- 
ble for distances of at least 5 A ]l3] , p2[ . Di Ventra and 
Pantelides employed density functional theory to deter- 
mine the relaxation of surface atoms for AhllO) surfaces, 
mimicking the tip by a single Al atom [|23fl . Within the 
range of distances from the point contact to the tunneling 
regime (1.3 to 3.3 A) they predict no corrugation due to 
atomic relaxations. The present calculation is in a sense 
complementary to their work. It is shown here that ex- 



perimental conditions in STM scans do not allow such a 
close approach, because the system becomes unstable al- 
ready at larger distances (4.6 A). Our calculations estab- 
lish an early onset of relaxations, which depend strongly 
on the position of the tip. For the on-top position of the 
STM tip we observe an onset of forces and relaxations 
at about 4.7 A, about 1.0 A higher than suggested by 
earlier simulations f|||. The surface atom reaches its 
furthest elongation from the surface at a tip-sample sep- 
aration of 4.5 A; it is then 1.3 A above its position on 
the isolated surface. This point marks the beginning of 
atomic transfer from the surface to the tip, the origin of 
the hysteresis observed in single approach/retraction cy- 
cles H . Neighboring atoms on the surface are equally re- 
laxed, their maximum displacement is 0.1 A. The second 
nearest neighbors remain close to their original position 
(vertical displacement less than 0.04 A). Even though the 
boundary conditions at our two dimensional repeat unit 
constrain the relaxations of individual atoms, the very 
small relaxation of second nearest neighbors indicates a 
rapid decay of the forces in lateral direction. Therefore 
we do not expect substantial changes in a calculation per- 
formed with a larger unit cell. If the tip is in the threefold 
hollow, the onset offerees and relaxations is substantially 
retarded and occurs at 4.3 A. The maximum outward re- 
laxation of the three gold atoms is somewhat lower and 
about 1.0 A (reached when the tip-sample separation is 
4.1 A). 

The tungsten tetrahedron of the STM tip relaxes about 
the same amount in both positions: the apex atom ap- 
proaches the sample by about 0.5 A at the point of clos- 
est approach. The change is due to an increase of the 
interlayer distance of the tungsten atoms by about 0.1 
to 0.2 A, and a buckling of the layers underneath, which 
extends about two layers into the lead. The forces on 
the tungsten tip are shown in Fig. [l]. The values for 
tip-sample relaxations are given in Fig. ^|. Owing to 
the strong interactions between tip and sample at close 
range, the core-core distance between tip and sample is 
decreased. The net effect of relaxations is a decrease 
of the distance between tip and sample. A deviation of 
about 2 A between the distance inferred from simula- 
tions, and the distance inferred from measurements of 
the decay, has actually been observed [Q. In the crit- 
ical range of 4.6 to 4.1 A the actual distance changes 
by about 2 A. The real distances are shown in Fig. || 
(a). This points to an effect frequently observed in ex- 
periments with atomic resolution: it is difficult to obtain 
atomically-resolved images, because sizable corrugations 
in most cases require a distance below 5 A, but it is very 
easy to destroy an STM tip by small increases of the 
current. We now understand that the actual range of 
imaging is less than 1.0 A, or less than an order of mag- 
nitude in the current scale. Stability is destroyed, when 
the apex atom of the tip and an atom of the sample form 
a chemical bond. To check, whether a different metal sur- 
face would lead to different results and thus a different 
interval in actual measurements we computed the forces 
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between a tungsten atom and a Cu, Ag, and Au atom 
in a dimer configuration. The result of our calculation 
is shown in the inset of Fig. [l| It can be seen that the 
forces are equal for all three noble metals in the range 
above 4 A. Short of performing identical simulations for 
all three surfaces this seems the best confirmation that 
our result is of general validity. 

Currents were calculated in a perturbation approach 
pl[ and based on the fully relaxed sample and tip sys- 
tems. The current was calculated for a sample bias of -0.1 
V. The curves in Fig. [3] (b) show the current for the hol- 
low position. We obtained the change due to surface and 
tip relaxations by calculating the full range of currents 
in a frozen atomic configuration. The logarithmic graphs 
are different for every distance, the z-labels give the tip- 
sample separation before relaxation. The true increase 
of the current is determined by taking only a single point 
of each graph. In practice we interpolated between the 
computed points to obtain the smooth black curve shown 
in the figure, giving an enhanced increase in the current. 
For an experimenter who estimates the potential barrier 
between tip and sample from the exponential decay of 
the current, the apparent barrier height increases in this 
range |l7],|25|]. The current calculation is based on the 
Baardeen approach p6fl , where sample and tip are con- 
sidered decoupled. In a previous paper we studied the 
effects of a tip potential on the surface electronic struc- 
ture and it was established that these effects are 
limited to very low distances (below 4 A). Here we find 
that the onset of relaxations begins about 1 A before 
this point. The main effect, not included in perturbation 
theory, is therefore due to relaxations. But this indi- 
cates that a more refined treatment of STM currents e.g. 
by a scattering approach is no substantial improvement 
compared to the Bardeen approximation, if it does not si- 
multaneously treat relaxations. And this, in turn, is still 
computationally untractable because it requires a fully 
ab initio treatment for every single position of the tip. 

All the results presented here are derived from ab ini- 
tio treatments of currents and relaxations for two lateral 
positions of the tip: the on-top position and the threefold 
hollow. The tunnel currents in our simulations were as- 
sumed to be 5.1 nA in both cases, but owing to the site- 
dependent relaxation of the systems the two constant- 
current contours differ. The piezoscale would give a dis- 
tance of 4.7 A on top of a gold atom and 4.5 A in the 
fee hollow site. These distances are at the lowest limit of 
stability for the tip-sample system. The difference of 0.2 
A describes the experimental corrugation of the Au(lll) 
surface, when measured with a tungsten terminated tip 
. The simulated scans with the tips and surfaces frozen 
in the two relaxed configurations are shown in Fig. |] (a), 
and (b). We stress that these images are idealizations, 
because they were computed with the wavefunctions of 
frozen atomic positions of tip and sample. In Fig. |^ (c) 
we show two simulated linescans across a single atom. 
The two separate curves describe the linescan for the two 
setups used (on-top, and hollow). It can be seen that the 



relaxation leads in effect to an enhancement of the corru- 
gation by about 0.13 A. Only with this relaxation do we 
obtain agreement between experiments and simulations 

The most remarkable result of our calculations, and the 
one that experimenters may find most important, is the 
very low range of stability for atomically resolved STM 
scans on metal surfaces. This range is only about 0.5 A. 
Below this point, at about 4.5 A from the surface, the 
distance between tip and sample changes drastically by 2 
A (within 0.5 A of measured displacement), as the two 
separate systems jump into contact. In this range the 
current may even decrease as the tip approaches, if one 
atom of the tip or sample system becomes separated. 
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FIG. 1. Forces on the apex atom (full triangles) and the 
subsurface atoms (empty triangles) of the tip above the hollow 
position of the Au(lll) surface. The forces on the apex atom 
are about one order of magnitude higher. The range of forces 
is about 4 A. Inset: forces on Cu, Ag, and Au atoms in a 
dimer. The forces between the tungsten atom and the noble 
metal atom are roughly equal in a distance range above 4 A. 
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FIG. 2. Relaxed positions of sample and tip atoms 
for three piezoscale readings of the STM: 5.0 A(left), 4.5 
A(center), 4.0 A(right); and two STM tip positions. At the 
point of closest approach the closest Au surface atoms are 
about 1 Ahigher than on the isolated surface. Note that re- 
laxations differ substantially for the two different positions. 
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FIG. 3. Relaxations and currents of the coupled system. 
The separation decreases, within a distance of 0.5 A by as 
much as 2.0 A (a). Logarithmic current graphs computed 
from the wavefunctions of the relaxed system. Individual 
curves pertain to equilibrium positions at a specific piezoscale 
reading. In the range of measurements at about 4.5 - 5.0 A 
the current decays faster due to relaxation effects, the actual 
barrier height seems increased (b). 
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FIG. 4. Current contour plots using the relaxed configu- 
rations of sample and tip for a distance of 4.7 A (tip on top of 
an Au atom, frame (a)), and 4.5 A (tip in the threefold hollow 
position, frame (b)). Simulated linescan across a single atom 
(c), for the setup in the on-top position (black curve, A- A' 
in frame (a)) and for the setup in the hollow position (grey 
curve, B-B' in frame (b)). All scans correspond to a current 
of 5.1 nA. An actual scan would have enhanced corrugation 
due to relaxation effects (see dash dotted lines in frame (c)). 
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